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ABSTRACT 

Weak gravitational lensing by the intervening large-scale structure of the Universe induces high-order 
correlations in the cosmic microwave background (CMB) temperature and polarization fields. We con- 
struct minimum variance estimators of the intervening mass distribution out of the six quadratic com- 
binations of the temperature and polarization fields. Polarization begins to assist in the reconstruction 
when £^-mode mapping becomes possible on degree-scale fields, i.e. for an experiment with a noise 
level of ~ 40/iK-arcmin and beam of ~ 7', similar to the Planck experiment; surpasses the temperature 
reconstruction at ~ 26/iK-arcmin and 4'; yet continues to improve the reconstruction until the lensing 
5-modes are mapped to / ~ 2000 at ~ 0.3/iK-arcmin and 3'. Ultimately, the correlation between the 
E and B modes can provide a high signal-to-noise mass map out to multipoles of L ~ 1000, extending 
the range of temperature-based estimators by nearly an order of magnitude. We outline four applica- 
tions of mass reconstruction: measurement of the linear power spectrum in projection to the cosmic 
variance limit out to L ~ 1000 (or wavenumbers 0.002 ^ /c ^ 0.2 in /i/Mpc), cross-correlation with 
cosmic shear surveys to probe the evolution of structure tomographically, cross-correlation of the mass 
and temperature maps to probe the dark energy, and the separation of lensing and gravitational wave 
5-modes. 

Subject headings: cosmic microwave background - dark matter — large scale structure of universe 
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Th j weak gravitational lensing of cosmic microwave 



^-ni ode p olarization (Kamionkowski, Kosowsky fc Steb _ 
bin^ 1997; [Zaldarriaga fc Seljak] 1997). Lensing converts 



background (CMB) temperature and polarization aniso- 
tropics provides a unique opportunity to map the distri- 
bution of matter on large scales and high redshift where 



densi W fiuct nations w^crc still linear. Although lensing cf - 



fects are app arent in the p ower spectra of tempe rature and 
polarization ( jSeljak 1996; Zaldarriaga fc Seljak 1998), it is 
the higher order correlations i nduced by le nsing that make 
mass reconstruction possible (Bernardeau 1997). 

By remapping the CMB fields according to potential 
gradients, lensing acts as a convolution in Fourier space 
which introduces correlations between angular wavenum- 
bers or multipole moments. From a quadratic combination 
of the multipoles, one can form estimators of the poten- 
tial fie ld and henc e the interve ning mass. [Zaldarriaga fc 
Seljak (1999) and puzik et al.| (2000) constructed noisy 
estimators out of the product of gradients of the tempera- 
ture and polarization fields. |H^ (2001a, b) showed that the 
minimum variance estimator constructed from the temper- 
ature field has substantially greater signal-to-noise with 
arcminute resolution CMB maps. This estimator enables 
mapping of the dark matter above the degree scale, where 
the defiection power peaks. The cosmic variance of the 
CMB temperature field itself prevents mapping on smaller 
scales. In this Paper, we show how this limitation can 
be overcome with minimum variance quadratic estimators 
involving the polarization field. 

The CMB polarization field in principle provides a more 
direct probe of lensing than the temperature field. Un- 
like the temperature anisotropy, there is negligible cos- 
mological contami nati on of the polarization field in the 
arcminute regime (Hu 2000a). Furthermore, density per- 
turbations in the linear regime generate only the so-called 



£^-mode p olarization to its comp lement, the 5-mode po- 
larization ( Zaldarriaga fc Seljak 1998). Although gravita- 
tional waves also generate 5-mode polarization , they do 
so on ly above the degree scale. [Hu| (2000b) and Benabed] 
et al.| (2001) used the induced correlation between the E 



and the 5-modes to construct statistical measures of the 
lensing. Here we show that the correlation allows a direct 
reconstruction of the lensing masses which in fact has in 
principle the highest signal-to-noise of all the quadratic 
estimators. 

In practice, achieving this potential in the presence of 
detector noise, systematics and foreground contamination 
of the polarization will be challenging. These same chal- 
lenges will also have to be overcome in order to probe the 
physics of the early universe through gravitational wave 
5-mode polarization. A lensing study can therefore be 
conducted as secondary science for an experiment devoted 
to gravitational waves. In fact, as the leading cosmological 
contaminant of the gravitational wave 5-modes, a lensing 
study may well be required of such an experiment 
2001c). 

We begin in §|| with a brief review of lensing effects on 
the temperature and polarization fields. In §||, we present 
a formal study of the minimum variance quadratic estima- 
tors of the lensing potential and show that the EB com- 
bination can produce a high signal-to-noise mass map out 
to the 10' scale. We explicitly construct this estimator in 
§^ and simulate its performance in the presence of detec- 
tor noise. We discuss applications of mass reconstruction 
in §11 For illustrative purposes, we use a fiat ACDM cos- 
mology throughout with parameters Vtc = 0.3, fib = 0.05, 
ftA = 0.65, h = 0.65, n = 1, Sh = 4.2 x 10"^ and no 
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Fig. 1. — An exaggerated example of the lensing effect on a 10° x 10° field. Top: (left-to-right) unlensed temperature field, unlensed 
^^-polarization field, spherically symmetric deflection field d(n). Bottom: (left-to-right) lensed temperature field, lensed ^^-polarization field, 
lensed ^-polarization field. The scale for the polarization and temperature fields differ by a factor of 10. 



gravitational waves. 



2. LENSING 



Weak lensing by the large-scale structure of the Universe 
remaps the primary temperature field 0(n) = AT(n)/T 
and dimensionless Stoke s para meters Q(n) and U (h) as 
(iBlanchard fc Schneider 1987; Bernardeau 1997; |Zaldar-| 
riaga fc Seljakl 1998) 



e(n) 

[Q±iU]{h) 



6(n + d(n)), 
[g±i/7](n + d(n)), 



(1) 



where n is the direction on the sky, tildes denote the un- 
lensed field, and d(n) is the defiection angle. It is related 
to the line of sight projection of the gravitational potential 
^(x,i:>) as d = V(^, 



(n) 



■I 



dD 



{Ds - D) 
DD. 



^{Dh,D), 



(2) 



where D is the comoving distance along the line of sight in 
the assumed fiat cosmology and Ds denotes the distance to 
the last-scattering surface. In the fiducial cosmology the 
rms defiection is 2.6' but its coherence is several degrees. 

We will work mainly in harmonic space and consider suf- 
ficiently small sections of the sky such that spherical har- 
monic moments of order (/,m) may be replaced by plane 
waves of wavevector 1. The all-sky generalization will be 



presented in a separate work (Okamoto & Hu, in prep). 
In this case, the temperature, polarization, and potential 
fields may be decomposed as 



e(n) 



dH 

(27r)2 

dH 



e(i)e 
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[Q±iU]in) = - J AL[Eil)±iBil)]e 
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where ipi = co s ■'■(x • 1). Lensing changes the Fourier mo- 
ments by P3 2000b) 



56(1) 



dH' 
(27r)2 
72;/ r 



e{i')w{i',L), 



(4) 



SE{1)= J ^^[^(l')cos2^,a-B(l')sin2^ia]w^(l',L), 



SBil) = 



dH' 



B{\') cos 2ipv\ + £;(!') sin 2(/?ia W{\', L) , 



(27r)2 

where i^yx = i^y — ip\,'L = \ — \' , and 

W^(1,L) = -[1-L]</.(L). 



(5) 



Here JO = O — O for example. In Fig. we show a toy 
example of the effect of le nsing on the tem perature and po- 
larization fields (see also Benabed et al.| 2001). The effect 
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Fig. 2. — Power spectra of the CMB temperature and polarization 
fields compared with the detector noise of the Planck satellite and 
a nearly perfect experiment with a noise level of = Ap/\/2 = 
1/iK-arcmin and a beam of cr = 4' (long dashed lines, thick for 
polarization, thin for temperature). The Planck experiment has 
sufficient signal-to-noise to map the B field but can only marginally 
map the ^^-polarization field; the nearly perfect experiment can map 
the all three fields to / = 2000. 



on the ^^-polarization is similar to that of the temperature 
and reflects the fact that cos2(p\'\ 1 for L <C where 
the lens is smooth compared with the field. Even in the 
absense of an unlensed 5-polarization, lensing will gener- 
ate it. The lensing structure differs since sm2(py\ ^ for 
L <^ 1. This fact will ultimately lead to a different range 
in L of sensitivity to (j) from the various fields. 

Since the unlensed fields and potential perturbations are 
assumed to be Gaussian and statistically isotropic, the sta- 
tistical properties of the lensed fields may be completely 
defined by the unlensed power spectra 



{x*mi')) 



{2-K fS{l - \')Cr' , 
(27r)25(L - L')L-^C 



where x = &, E, B and we have chosen to express the 
potential power spectrum with a weighting appropriate for 
the deflection field d{ii). Under the assumption of parity 
invariance 



C, 



EB 



0, 



(6) 



and in the absence of gravitational waves and vorticity 
= 0. The peak in the logarithmic power spectrum 
L2Cf^/27r at L - 30-40 defines the de gree-scale coherence 
of the deflection angles. 

Finally, we define the power spectra of the observed tem- 
perature and polarization fields as 



{x*{V)x(l'))^{27rfS{\-l')C[ 



(7) 



where the power spectra include all sources of variance to 
the fields including detector noise and residual foreground 
contamination added in quadrature. We will inc lude Gaus- 
sian random detector noise of the form (Knox 1995) 



\ ^l{l^l)a^/8ln2 

^CMB 



EE 



BB 



\ TcMB ) 



^^(Z+l)cr2/81n2 



(8) 



where A^^p parameterizes white detector noise, here in 
units of /iK-radian, Tcmb = 2.728 x 10^ /iK, and a is the 
FWHM of the beam. We will often assume Ap = V^At 
as appropriate for fully-polarized detectors. In Fig. |[ we 
compare the signal and noise contributions to the total 
power spectra for the Planck satellite experiment^ (mini- 
mum variance channel weighting from Cooray & Hu 2000; 
At ~ 27/iK-arcmin, Ap ^ 40>/2 /iK-arcmin, a ^ 7 ) and 
a near perfect reference experiment (At = Ap/V^ = 
1/iK-arcmin and a = 4'). In general where the signal 
exceeds the noise power spectrum of a field, there is suf- 
ficient signal-to-noise for mapping. When this is not the 
case, a statistical detection of the signal may still be pos- 
sible. The Planck experiment is on the threshold of being 
able to map the ^^-polarization. The reference experiment 
can map all 3 fields to / ~ 2000. 

3. MINIMUM VARIANCE ESTIMATORS 

As can be seen from Eqn. (H), lensing mixes and there- 
fore correlates the Fourier modes across a range defined by 
the power in the deflection field Cf^ (Hu 2000b). Consider 
averaging over an ensemble of realizations of the temper- 
ature and polarization fields but with a fixed lensing field. 
The two-point correlation of the modes takes the form 



(a;(iy(l'))cMB = /a(l,l')'/'(L), 



(9) 



where x, = 6, and L = 1 + 1'. We have assumed 

1^—1' and will use the subscript a to distinguish be- 
tween choices of the xx' pairing, e.g. a = 00. The cor- 
relation returns the value of the deflection potential with 
weightings fa that depend on the unlensed power spectra 
of Eqn. (M), which are given explicitly in Tab. 1. 

The two point correlations of the CMB Fourier modes 
themselves cannot be used to reconstruct the deflection po- 
tential since (j) is also statistically isotropic so that in the 
true ensemble average {(1>{L)) = 0. Eqn. does suggest 
however that an appropriate average over pairs of multi- 
pole moments can be used to estimate the deflection field 
d(n). 

^http: / /astro. estec.esa.nl/Planck 



a /a(ll,l2) 



ee 
eB 

EE 
EB 
BB 



C'Q^cos^,,,,(L.li) + C'Q^(L.l2) 

C'«^sin2^,,,,(L-li) 

[C;f^(L-li) + C';f^(L-l2)]cos2^,,,, 

[C,f^(L.10-C',f^(L-l2)]sin2^.,,, 

[C',fS(L-li)+C,f^(L-l2)]cos2^n,, 

Table 1 
Minimum variance filters 
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Let us define a general weighting of tlie moments 



AaiL) f d% 



L J (27r)2 
where I2 = L — li and the normalization 



x{h)x'{h)F^{h,h), (10) 



Aa{L) = L' 



is chosen so that 

{da(L))cuB = d{-L) = L4>(L) . 



(11) 



(12) 



In general there are 6 estimators corresponding to the 3! 
pairs of B. In the assumed cosmology, where gravi- 

tational wave perturbations are negligible compared with 
density perturbations, a = BB has vanishing signal-to- 
noise effectively reducing the estimators to 5. 

We can optimize the filter Fq^ by minimizing the variance 
{dl^{L)da{L))^ subject to the normalization constraint 



Fa{h,h) 



^xx ^x' x' ^x'x' ^xx 
h h h h 



^iQxx' ^xx' ^2 



(13) 

This filter takes on simple forms for two common cases: if 
X = x' , as in the case of a = QQ, EE and BB, 



Faih,h) 



2(jxx^xx ' 



(14) 



if q 



; — 0, as in the case of a = OB and EB, 



^XX ^X X 

The noise properties of these estimators follows from 



(15) 



(d;(L)d^(L')) = (27r)25(L - L')[Cf + iV„^(L)] , (16) 



where 



FMl2,li)Cr;"^cf^''). (17) 



Recall that the xx-power spectra account for both the cos- 
mic variance of the fields and the noise variance of the 
experiment. Notice that for the minimum variance filter 



NUL)=AaiL). 



(18) 



In Fig. |[ we compare the signal and noise power spectra 
for the Planck experiment and the reference experiment 
defined in §^ Recall that true mapping is possible when 
the signal exceeds the noise spectrum. For the Planck 
experiment, 00 provides the best estimator reflecting the 
fact that Planck will not be able to produce true maps of 
the polarization modes. Furthermore, the signal-to- noise 
is highest at L ^ 200 reflecting the fact the modes are 
mainly correlated across AL ~ 60, where the deflection 
power spectrum peaks. 




Fig. 3. — Deflection signal (dd) and noise power spectra of the 
quadratic estimators and their minimum variance (mv) combina- 
tion: (a) Planck experiment (b) reference experiment. As the sen- 
sitivity of the experiment improves the best quadratic estimator 
switches from BB to EB. Only the £^5-estimator can reconstruct 
the mass distribution at L > 200. 



For the reference experiment, all 5 estimators have suf- 
ficient signal-to- noise to produce maps at L < 200. The 
EB estimator has the best signal-to-noise, and allows for 
mapping to L ^ 1000. The reason is that there is no noise 
variance contributed by an unlensed B field. Furthermore, 
the signal intrinsically comes from higher L. A 5-field 
at a wavenumber 1 cannot be generated from neighboring 
modes 1' ~ 1 from the low L deflection field because of the 
sin term in the lensing kernel (see Eqn. ||). Thus the signal 
to noise is relatively higher at high L in the EB estimator. 

For experiments that are intermediate in sensitivity be- 
tween Planck and the reference experiment, the five es- 
timators of the deflection field have comparable signal- 
to-noise and may be used to cross check each other. At 
high-L where the individual estimators are noise limited, 
combining the estimators as 



G^mv(L) 



Va{L)da{'L) , 



(19) 



can substantially reduce the noise. The minimum vari- 
ance weighting is a generalization of the inverse variance 
weighting that accounts for the covariance in Eqn. dlj 



E/3^(N-l)/5, 



(20) 



Fig. 5. — Mass reconstruction on a 10° x 10° field with the reference experiment (At = Ap/\/2 = 1/iK-arcmin and a = A'): (a) defiection 
field, (b) BO-reconstruction, (c) £^B-reconstruction. 




Fig. 4. — Defiection signal {dd) and noise power spectra for 
the minimum variance (mv; solid lines) and BB (dashed lines) 
estimators as a function of (a) beam size cr and (b) noise level 
At = Ap/ \/2. The noise saturates to its minimum by cr ~ 2 — 4^ 
and At ~ 0.1 — 0.3/iK-arcmin as the polarization field is mapped 
to the cosmic variance limit out to / < 2000. 



The noise variance 

(Cv(L)t;mv(L')) = (27r)25(L - L')[Cf + N^,{L)\ (21) 
becomes 



Note that the 00 and EB estimators are independent and 
those estimators that are correlated have correlation coef- 
ficients Ncc^l ^^NoccxNl3j3 of no more than tens of per cents. 

The minimum variance noise spectra for Planck and the 
reference experiment are shown in Fig. I We give it as 
a function of the noise At and beam a in Fig. ||. The 
signal-to-noise saturates around At ^ 0.1 — 0.3/iK-arcmin 
and (J ^2- A'. Below At ^ 26/iK-arcmin {/\p ^ 37/iK- 
arcmin), the combined signal-to- noise in the QE and EE 
estimators exceeds that in 00 at L 40 and a = A' . For 
the smaller scales oiL ^ 300, where only the EB estimator 
plays a role, the EB and 00 estimators have comparable 
signal-to-noise around At ~ 6/iK-arcmin (Ap 8/iK- 
arcmin) at L ^ 300 and <j = 4'. 

4. EB ESTIMATOR 

As we have seen in the previous section, the EB- 
estimator of the defiection field has the potential to map 
the mass distribution out to L ^ 1000. We therefore ex- 
plicitly construct and test this estimator in this section. 
This construction i s ve ry similar to that of the 00 esti- 
mator presented in Hu (2001b). 

From Eqn. (p^Oj), the EB estimator is 

(23) 

where recall L = 1 + 1'. The convoluted form of this esti- 
mator suggests that it may be re-expressed as a product 
of fields on the sky. To see this, rewrite 

sin2(/?ii/ = 2(l-r)[n- (1 x 1')], (24) 
where n = —63. We can then define the filtered fields 

E^jk{-^)= / ^^l{lihh)^E{\)e''-'^ , (25) 



(27r)2 

dH 



1 



ilfi 



(27r)2^"'"^'C,«S' 



(26) 



E„.(N-^)a/3 



(22) 



There are 4 unique filtered £^-fields and 3 unique filtered 
5-fields. They may be combined to form the appropriate 
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FlG. 6. — Statistical errors achievable on the deflection power 
spectrum with the Planck (/sky = 0.65) and reference experiments 
(/sky = !)• Boxes represent band averaging width and la errors. 
The polarization information in the reference experiment allows for 
a cosmic variance limited measurement of the projected power spec- 
trum out to L ~ 1000. In this regime, the fluctuations are almost 
completely linear (dashed lines). 



dot and cross products 



Gi{h) = 2 ^ Eijk{n)Bjm{n)ekm3 

jkm 



(27) 



where Cijk is the Levi-Civita symbol. The deflection field 
is then reconstructed as 



L-G(L) 



(28) 



The other quadratic estimators can be constructed in a 
similar fashion. 

Fig. H shows an example of the EB reconstruction 
compared with the 00 reconstruction on a 10° x 10° 
field with the reference experiment. Notice that the EB- 
reconstruction has substantially lower noise on small angu- 
lar scales. We assume here that the unlensed power spec- 
tra have been determined externally from precision satel- 
lite missions and through the modelling with cosmological 
parameters (see Hu 2001b). Errors in the determination 
translate into non-optimal filters and a small bias in the 
amplitude of the reconstructed maps. 

5. APPLICATIONS 

In this section, we outline four applications for mass re- 
construction: measurement of the (linear) power spectrum 
in projection, cross correlation with cosmic shear observa- 
tions, cross correlation with the temperature field, and 
decontamination of the polarization signature of gravita- 
tional waves. The fir st t hree applications have been ex- 
tensively discussed in Hu (2001c) for the 00 temperature 
based estimator and we refer the reader to details therein. 
Here we focus on the additional information provided by 
the polarization field. 

5.1. Linear Power Spectrum 

The most direct application of mass reconstruction is 
to measure the matter power spectrum in projection, i.e. 
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Fig. 7. — Statistical errors on the cross correlation of CMB deflec- 
tions and cosmic shear in three source redshift bands on a 1000 deg^ 
patch of sky for the Planck and reference experiment. Assumptions 
for the cosmic shear experiment are given in the text. Precision 
measurements from the polarization estimators enables highly sig- 
niflcant cross correlation detection and hence tomographic studies 
of structure evolution. 



the deflection power spectrum Cf^^ itself. Power spectrum 
measurement requires only a statistical detection of the 
deflection field, not a true reconstructed map and there- 
fore can be extended to higher wavenumbers or smaller 
scales than is possible for mapping. The noise level for the 
estimation of band powers is reduced by averaging over L 
directions in a band AL 



LAL/sky 



iVn,v(i)]. 



(29) 



where /sky is the fraction of the sky covered by the ex- 
periment. In this approximation, the noise is assumed 
to be Gaussian. This should be a good approximation 
where the sample variance of the lenses dominates the 
noise variance. Formally, the noise will be increasingly 
non- Gaussian at high L as the estimator is constructed 
out of fewer arcminute scale temperature and polarization 
fluctuations. Quantification of this effect for the temper- 
ature based reconstruction show that its effects are minor 
( |Hu| 2001a); a full treatment requires the consideration 
of the temperature-polarization trispectrum (Okamoto & 
Hu, in prep.). 

Polarization enables two advances over what can be 
achieved by the temperature field alone. As in the case 
of mapping, polarization enables precision measurements 
at small scales through the £^5-estimator. In Fig. 
we compare the Planck experiment (with /sky = 0.65 
and the reference experiment (with /sky = 1); as seen in 
Fig. ^, former relies mainly on the 00-estimator and the 
later on the £^5-estimator. The noise in the Gaussian 
approximation approaches the sample variance limit of 
ACjf/Cjf = (LAL/sky)~^/^ on the scales L ^ 1000, i.e. 
a total of 1% precision in each 1% of sky. This corresponds 
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Fig. 8. — Statistical errors on the cross correlation of CMB de- 
flections and the temperature fleld for the Planck (/sky = 0.65) and 
reference (/sky = 1) experiments. The flve estimators of the de- 
flection fleld obtainable with polarization information enables flve 
nearly independent, cosmic variance limited detections of the cross 
correlation for L < 100 (shown offset slightly for clarity). The cross 
correlation is extremely sensitive to the properties of the dark en- 
ergy. 




to scales in the matter power spectrum of 0.002 < /c < 0.2 
in /i/Mpc representing the whole linear regime today. 

Equally importantly, the polarization allows for sharp 
consistency tests on the power spectrum measurements at 
L <^ 500. In the reference experiment, all 5 estimators 
have sufficient signal-to-noise to measure the power spec- 
trum here. It is highly unlikely that any unknown contam- 
inant from foregrounds or instrumental systematics would 
affect specific quadratic combinations of the temperature, 
^^-polarization, 5-polarization, in the same way. 

5.2. Evolution of Structure and Cosmic Shear 

One would like to go beyond the projected power spec- 
trum to the three-dimensional distribution to track the 
evolution of structure and hence the physical properties of 
the dark matter and energy. This is not possible through 
CMB lensing alone since the source plane lies at the ef- 
fectively infinite redshift of last-scattering. Weak lensing 
also distorts the shape of distant, but for our purposes 
foreground, galaxies allowing a measurement of the gradi- 
ent of the defiection angles, or more properly, the s o-called 
cosmi c shear, from wide-field ini aging surveys (see Mellier 
1999; Bartelmann & Schneider 2001 for reviews). Since 



these sources are distributed across a range of redshifts, 
the change in the mass reconstruction as a function of 
source redshift can p robe the radial d istribution of mat- 
ter tomographically ( Wittman et al. 2001). The effect 



on the power spectra and cross-correlation of cosmic shear 
C|f has been shown to be an effec tive probe of the dark 



energy equation of state ( Huterer 2001; Hu 2001c). Be- 
cause cosmic shear studies are most effective for L ^ 100, 
tomography with the lensing of the CMB temperature is 
difficult. 

By extending the measurements to overlapping wavenum- 
bers, CMB polarization allows tomographic studies to be 
anchored at the high-z end. In Fig. 0, we show the errors 
on the CMB defiection-cosmic shear cross power spectrum 
achievable with the Planck vs. reference experiment 




Fig. 9. — Large-angle (/ < 100) lensing B-polarization fleld (top) 
and the reconstructed B-polarization fleld from the small angle EB 
deflection estimator and the observed £J-fleld. Detector noise appro- 
priate for the reference experiment has been added to this 25° x 25° 
patch. Reconstruction techniques can help separate the gravita- 
tional wave and lensing S-modes. 



and 1000 deg^ of overlap with a cosmic shear survey out to 
median source redshift z = 1^ divided into three redshift 
bands < 1, 1 < < 1.5, > 1.5. Errors on the cosmic 
shear side assume n = 56 gal/arcmin^, and an intrinsic 
shear measurement error of (7int)^ = 0-4 per component 
per galaxy. 

5.3. Dark Energy and the ISW Effect 

The integrated Sachs- Wolfe (ISW) effect from the differ- 
ential redshift due to the decay in the gravitational poten- 
tial is extremely sensitive to t he bac kground proper ties of 
the dark energy ( [Coble et ah] 1997; [Caldweh etal] 1998) 
and provides a unique handle on its clustering properties 
( |llu| 1998). The latter can potentially test the scalar-field 
hypothesis for its nature. Unfortunately the ISW effect is 
buried under the larger primary anisotropy and can best 
be isolated through cross correlation with other large-scale 
tracers of the gravitational potential. The defiection field 
of CMB l ensing provides a perf ect ca ndidate for cross cor- 
relation ( Seljak fc Zaldarriaga] 1998, Goldberg fc Spergei] 
1999, Hu 2001c). In a flat universe, detection of any cross- 
correlation at all represents an essentially direct detection 
of the dark energy. 
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Since the cross correlation effect is confined to multi- 
poles L <^ 100, the 00 estimator can itself reach the cos- 
mic variance limit for detection. What polarization pro- 
vides is 4 other nearly independent probes of the cross- 
correlation. Furthermore the polarization estimators each 
contain enough signal-to-noise to reconstruct the deflec- 
tion field with independent sets of multipoles in the polar- 
ization field (see [Hu| 2001 for the analogous technique using 
the temperature field). Since the signal is weak and the 
importance of understanding the particle properties of the 
dark energy great, the ability to make these consistency 
tests is an important asset. 

In Fig. ^, we compare the ability of the Planck experi- 
ment (/sky = 0.65) and the reference experiment (/sky = 1) 
to measure the deflection-ISW cross correlation Cf^. 

5.4. Gravitational Waves and B-modes 

The 5-mode polarization produced by gravitational 
waves offers what i s perhaps our most direct window on 
the ea rly universe (p^amionkowski, Kosowsky fc Stebbine 



1997; ^aldarriaga fc Seljak 1997). Under the inflation 
ary paradigm, its amplitude determines the energy scale 
of inflation. In this context, gravitational lensing, which 
also generates 5- modes acts as a contaminant. To con- 
strain inflationary energy scales below 10^^ GeV, removal 
of the lensing contaminant, either statistical or direct, will 
be required (Hu 2001c). Fortunately the converse prob- 
lem does apply: since the 5-modes used to reconstruct 
the deflection fields reside in the arcminute regime, even 
much larger amplitude gravitational wave 5-modes at de- 
gree scales do not contaminate mass reconstruction. This 
fact allows the possibility of direct removal of the lensing 
5-modes at degree scales. 

We save detailed exploration of 5-mode decontamina- 
tion to a future work and here simply show that the mass 
reconstruction from the small-scale 5-polarization itself 
has sufficient signal-to-noise to make the procedure feasi- 
ble. Decontamination is not feasible through the 00 es- 
timator since most of the 5-modes on degree scales arise 
from fluctuations in the deflection field with L ^ 200. 

Consider the £^-modes of the observed, i.e. lensed and 
including detector noise, polarization field and construct 
the Stokes parameters Q(p) and U{u) from it alone via 
Eqn. (H). Use the reconstructed deflection field dEsi'^) to 
artificially lens the distribution. The 5-field of the result- 
ing polarization is an estimator of the 5-field from lensing 
that is independent of the true 5-field on large scales. In 
Fig. |[ we show the true 5-field from lensing, low pass fil- 
tered to / < 100, compared with the reconstructed 5-field 
for the reference experiment on a 25° x 25° field. 

6. DISCUSSION 

Based on the induced correlation between the E and B 
modes, the lensing of CMB polarization offers the oppor- 
tunity to reconstruct the mass distribution in projection 
on scales corresponding to 0.002 < k < 0.2 in h Mpc~^. 
Compared with a similar reconstruction from the temper- 
ature field, polarization allows for an order of magnitude 
extension to smaller scales. These small scales can corre- 
spondingly be probed with the smaller degree scale fields of 
view that are more typical for planned polarization stud- 
ies. Moreover mass reconstruction is only sensitive to the 
arcminute scale correlations in the polarization field and 



does not require a true map over the full field. This ad- 
ditional range does not come at the expense of higher res- 
olution requirements: the signal-to-noise saturates at the 
several arcminute scale corresponding to / ~ 2000. Due to 
the smaller absolute scale of the signal, polarization stud- 
ies do require much more sensitive detectors, with the EB 
estimator surpassing the 00 estimator for Ap ^ 8/iK- 
arcmin for L ~ 300. Sensitivity and control over fore- 
grounds and systematics are issues that any polarization 
based study must address, especially those searching for 
the gravitational wave imprint from inflation. 

Mass reconstruction from CMB polarization can pro- 
vide measurements of the matter power spectrum over a 
wide range of scales that are entirely free of assumptions 
of how the luminous matter traces the mass (or bias) and 
the distribution of lensing sources, as well as largely free of 
non-linear corrections. It complements cosmic shear stud- 
ies by providing the deepest two-dimensional mass maps 
possible to anchor tomographic studies of the evolution of 
structure. It extends the shear-based lensing studies to 
near horizon-sized structures and therefore provides the 
opportunity to study the dark energy in its cross correla- 
tion with the ISW effect in the temperature field. Finally, 
since reconstruction only requires information from fine- 
scale correlations in the polarization field, it may be used 
to remove the lensing 5-modes on large-scales from any 
potential gravitational wave signal. These potential scien- 
tific returns may help justify the great experimental effort 
that will be required to map the CMB polarization field. 
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